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We calculate the next-to-leading order (NLO) as-corrections to the contribntions of the conden¬ 
sates (aG^) and (qq)^ in the cnrrent-current correlator of the hybrid current gq{x)'yuiFlf,^T°'q{x) 
using the external field method in Feynman gauge. After incorporating these NLO contributions 
into the Laplace sum-rules, the mass of the ^ light hybrid meson is recalculated using the 

QCD sum rule approach. We find that the sum rules exhibit enhanced stability when the NLO 
Qs-corrections are included in the sum rule analysis, resulting in a 1 ^ light hybrid meson mass of 

approximately 1.6 GeV. 


I. INTRODUCTION 


Soon after Quantum Chromodynamics (QCD) became established as the theory of strong interactions, the search for 
mesons with exotic quantum numbers was initiated. Important experimental progress has occurred in identifying 
potential candidates for exotic mesons. Among these are the two = 1 '' isovector states /5(1405) and /5(1600) 
which have been identified by two collaborations in the channels pir, rjir and rj'?: i I I- This gives the theorists 
an opportunity to compare their results for the exotic light hybrid mesons with the experimental observations. For 
instance, the widely cited mass prediction for the = 1 ^ state from lattice QCD lies around 1.9 GeV [Q which 
disagrees with the experimental data. The decays of the 1 ^ hybrid meson have also been studied in the context of 
various models an. and also appear to be in disagreement with the experimental data. A possible reason for these 
inconsistencies may derive from the fact that non-perturbative effects cannot be easily controlled at this low energy 
scale. A further possibility is, of course, that the two new states may have been misidentified and are not hybrid 
mesons after all. It is clear that further theoretical studies of the properties of the hybrid mesons are necessary before 
one can be confident of their discovery. 

In this paper, we concentrate on the mass prediction for the 1 ^ hybrid meson using the QCD Laplace sum rule 
approach [Q. Previous theoretical studies found masses in the range 1.4-2.1 GeV [^, with more recent estimates 
around 1.6 GeV [^, |^. These latter estimates are close to the /5(1600) but do not accommodate the state /5(1405) 
as a hybrid candidate. However, the QCD sum rule analysis is not very stable, which from sum-rule studies of scalar 
gluonium and the p meson, could be attributed to the importance of NLO corrections associated with the dimension- 
four (agG^) and dimension-six {qq)"^ operators |11, In particular, the four-quark operator has the same dimension 
as the two point correlator of the 1 ^ current and its coefficient function is non-logarithmic in the leading order. 
Thus the dimension-six contributions are absent in the next-to-lowest moment sum-rule, resulting in discrepancies 
with the mass prediction from the zeroth moment sum rule j^. Thus it is interesting to determine whether the 
Qfs-corrections can reduce this discrepancy. For consistency, the ag-corrections to both the dimension-four (agG^) 
and the dimension-six {qq)"^ operators need to be included. 

In this paper, we first give a brief introduction to the external field method used for our calculations. We then 
present our calculations for the Qfg corrections to {qq)'^ and {agC^). Finally, we discuss the effect of the ag-corrections 
on the mass prediction for the 1 hybrid meson. 


II. EXTERNAL FIELD METHOD 

In order to obtain the contribution of the operator asG“^‘^G“,^ in the Operator Product Expansion (OPE) expansion 
of the relevant current-current correlator, it is convenient to use the so-called external field method introduced in Ref. 
[ p^ . Because we will use a slightly modified version of the external field method, we give a brief introduction to this 
method before carrying out our calculations. We first split the gauge field into two parts 

= i^n)ext + ( 1 ) 

where is an external (classical) field and a“ is a quantum field. In the following we will suppress the label 

ext in (H^)ga,^ for convenience. Substituting ( 0 ) in the QCD Lagrangian one obtains 


( 2 ) 
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where 


^flL' — . [-^/i; 5 — dfl “1“ “t“ 

*5 


The Lagrangian (i is invariant under the transformation 


( 3 ) 


afj,{x)a'^ = U ^{x){af^{x) + A^{x))U{x) - Af_^{x) + ^{x)U{x), (4) 

ipix) ij)' = U{x)il){x), (5) 

Af, ^ A'^{x) = Af,{x) . (6) 

Unlike in |p^ , we demand that the external field is invariant under the gauge transformation. Therefore, any gauge 
condition imposed on the external field does not break the invariance under the gauge transformation This is 
consistent with the result given in |^, which states that the series of the OPE for a correlator of gauge invariant 
currents does not depend on the choice the gauge of the external field. Furthermore, in order to fix the gauge of the 
quantum field, we use the Feynman gauge 


-\{d^an^ (7) 

instead of the background gauge used in Q, where = dfj, + igA^. This has the advantage that 

it can be shown that the external method in the covariant gauge is equivalent to the plane-wave method [Q. Note 
that, when the external field method is used, the background gauge very likely differs from the covariant gauge in 
processes where both the radiation field and the external field are present in the initial states and (or) the final states. 
In the processes where only the radiation field or the external field is present in the initial states and (or) the final 
states, there is no difference between the background gauge and the covariant gauge. For instance, if there is only the 
external field present in the initial and (or) final states, the radiation field is integrated out. Then the external method 
in the background gauge is exactly the same as the original background gauge method ||^. The later is consistent 
with the plane-wave method. If there is only the radiation field present in the initial state and (or) final state, we can 
switch off the external field. Then there is no difference between the background gauge and the covariant gauge. 

The Feynman rules in the Feynman gauge can be obtained quite straightforwardly. Under the infinitesimal gauge 
transformation 

(5a“ + (8) 

one can easily derive the Lagrangian for the ghost field 

LgHost = - ^r\al + Al)]6,. (9) 

Then, the external field A^ obeys the same Feynman rules as the radiation field a^. This is also consistent with the 
plane-wave method |p^ . The calculational techniques using the Lagrangian (i are similar to those in the background 
gauge 

In order to extract the operator asG °'^'^from the OPE, the Fock-Schwinger condition is imposed on the external 
field 


x^A'^^ix) = 0 . ( 10 ) 

Then, with the aid of the technique proposed in [T^ , the gluon propagator in the presence of the external field can 
be obtained straightforwardly. For instance, up to the term 0{q~'^) the gluon propagator is given by 


D^iuiq,y)= / e*‘?'"^d^a; (0|ra“(a;)a()(y)|0) 


—i 


3G 


ah 


2' 


qP{G 


ah , 


G ab 
py 


9m) 


iy'' {Gpiqiy 


G' 




( 11 ) 


where and D = 4 — 2e is the dimension of space-time. As expected, the gluon propagator ( 0 ) differs 

from the corresponding one given in [^ . Sometimes it is more convenient to do the calculation in coordinate space. 
The gluon propagator in coordinate space can be obtained by using the Z?-dimensional Fourier transformation 




(27r) 


D ^ 


'iq,y) 


Df,u{x,y) 


( 12 ) 
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The necessary integration techniques in £)-dimension were given in |T^ . One can convert the 4-dimensional quark 
propagator given in to the _D-dimensional quark propagator 


r(f) ^ r(f -1) 

27r'? Stt^' 


(13) 


III. NLO 0{as) CORRECTIONS TO (G^) AND {qqf 
The two point correlator of the 1 ^ hybrid current is defined as 

= * J d4xe*^“(0|T{j—(x),jr+( 0 )}| 0 ) 

~ iQiJ.Qi' ~ 9fj,i^Q )n^((5 ) + )) 

where the current j^{x) = q{x)T°‘^yigF^jjq{x) carries isospin 7 = 1 and the invariants 11^(Q^) and ns((5^) correspond 
to the contributions from 1 and 0^“*' states respectively. The renormalized current is denoted by which in 

the massless quark limit, has the form 


fr = zj,. (15) 

where up to the order 0{as) and in the Feynman gauge the renormalization constant Z is given by 0 

2 o2 1 

The leading order calculation of (|T^ including the quark and gluon condensate contributions are contained in 
and the NLO corrections to the perturbative part of (0) were calculated in i 0- Next we consider the NLO 
corrections to the gluonic condensate {asG^^G^^'") contributions. We divide the calculations into two parts. One part 
can be obtained via the calculations of Feynman diagrams shown in Figure where we do not display diagrams which 
vanish in dimensional regularization or which can be obtained from diagrams a-o in Fig. by symmetry arguments. 
Because the expansion in term of violates translation invariance, the Fig. diagrams g, j, and k respectively differ 
from diagrams m, n, and o. A straightforward calculation results in the structure 

asGo.pG‘^^Iif {q^) + (^ - asGppGU^q^Iif {q^) 

+ [asG^gGt - ^asGppG^^qPq^^ nf (g") + asG^pG^^q^q'^uf (q^) 

+ (g2) + ^a.Gp^G^g^nf (g^). 

In order to extract {asGp^G^'^), we need a condition for the gluonic vacuum expectation value which reads 


(17) 


Gp.Gp^ - 


^a(3^ ^{.9fip9i'<y 9pio'9yp)' 


Then, in the MS-scheme and Feynman gauge, we obtain 


nf“(g^) = -:Tin(^)(a,G^) 


3071 




247r 


16 as(/r) 




sG^) 


9 TT 

16 as{n) 


In 


In 




1139 as{n) 


216 TT 
1459 as(/i) 


216 TT 


-|- infinite terms 
infinite terms , 


( 18 ) 

( 19 ) 

( 20 ) 


where {aG^) = Another part of the next-to-leading order calculation of {asG^iyG^^) results from the 


renormalization of the current (^), i.e. 

* J d\F‘>^2iZ-l){0\T{jp{x),j+{0)m, 


(21) 
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It reads 


n: 


G^b 


= Hr 


3071 


n 


G-'b 




) (H) 

(asG^) 


Qtt 


Qtt 



88 asifJ.) 
^27 TT 
104 asjfJ.) 
27 TT 


+ infinite terms 
+ infinite terms, 


( 22 ) 

(23) 


where the sum of the infinite terms in ([l^ ) and (|^ ) is scale-independent. We did not check on the infrared convergence 
of the sum of infinite terms because we used dimensional regularization. However, the sum of these two parts must be 
IR convergent so that the Wilson coefficient of the condensates only depends on short-distance effects. Obviously, this 
result is invariant if we use the background gauge, because all radiation fields are integrated out. We have checked 
such invariance. 



FIG. 1: Feynman diagrams for the Os-corrections to the condensate (osG^). Dots stand for the cnrrent vertices. 


Similarly, by calculating the Feynman diagrams shown in Fig. 0, we obtain the next-to-leading order corrections 
for the {qq)'^ contributions 


47r 


n?(g^) = j^asifi){qqf 

zItt 

= -^as{n){qqy 


11 , 1 \ asifJ-) . / -r 

+ F^/) 


72 6 

53 1 


72 ^ 6 ^^/ TT 


TT 


91 


108 18^^/ TT 




In 


-q 


14 1 


asifJ.) 


where we have used the vacuum saturation approximation 

/ ^b—a'^b'\ _ 7 1 \2/-\2 


{(liH' Vj') = {T7:) {qQ) 3'^abSa'b' “ Sij'k'jSab'5a'b) 


(24) 

(25) 

(26) 


IV. MASS OF THE 1’+ HYBRID MESONS 


QCD sum rules are based on the resonance plus continuum duality ansatz 

—/TOn„_s(s) = ^ M^f'^6 (s — M^) -I- QCD continuum, (27) 

^ R 

where M/j is the mass of the resonance R, fa denotes the coupling of the resonance to the current and we use the 
narrow resonance approximation (5(s — A7^) for the resonance R. The spectral density Pv,s{s) = ^ImJly^s{s) can be 
related to the correlator n^^s(( 7 ^) at the scale —q^ via the standard dispersion relation 


n.i;,s(g^) 





Pv,s 

:"(s - g2) 


n —1 




(28) 
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FIG. 2: Feynman diagrams for the Os-corrections to the condensate( 55 )^. Dots stand for the current vertices. 


where the Uk are appropriate subtraction constants to render Eq. (^ ) finite. The energy variable —q^ has to be chosen 
in a region where one can incorporate the asymptotic freedom property of QCD via the operator product expansion. 
The spectral function can be expressed as 




dim n 


(29) 


where the scale fi separates the long-distance and short distance regime of QCD. The Wilson coefiicients Cn{—q^,fJ,) 
for the low dimension operators were given in |^, ^ as summarized in Section 3 of this paper. Ref. also introduced 
a dimension two operator resulting from the resummation of the large order terms of the OPE series. However, this 
method appears to be model-dependent and could result in a double counting of the contribution of the operators 
considered in the present approach. We will therefore not include the dimension two term in our analysis. 

Concentrating on the analysis of the vector 1 ^ channel, the lowest-lying resonance in the spectral density is 
enhanced by the standard approach of applying the Borel transform operator to weighted by powers of [Q. 
This results in the Laplace sum-rules 


/ -So 


k = 0,1,2, 


(30) 


where the quantity Rk represents the QCD prediction, and the threshold sq separates the contribution from higher 
excited states and the QCD continuum. In the single narrow resonance scenario, the lowest-lying resonance mass can 
be obtained from ratios 


m 2 Rk+i ir, so) 

Rk {t, So) 


(31) 


The zero-weighted sum-rule Rq for n/ = 3 can be obtained from Eqs. (|19| ) and (P2D, and Refs.||§, |^, with some 
results needed for calculation of the necessary Borel transforms extracted from |]o. 


1 a 


1301a 


2407r2 TT V 240 tt 


1 a /17( 


i?o(DSo)=7rT^-(l + ^-) I / fV^Hn ( ^ ) dt 


1207r2 TT V 72 TT 


0 


0 


3071 V 72 TT 


i (i -1#?) + i (y;) I -i* 

0 0 


367r \ 9 TT 
Att / 47 a 


sq 

■y2\ I ^-tr 1 


(32) 


- (Oe) - Y (l + “(99)' ■ Y ( 


The quantity (Og) is defined by 




(33) 
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Renormalization-group improvement of is achieved by setting y? = \ and higher-weight sum-rules can be 

obtained from r derivatives of (|^ ) before implementing renormalization-group improvement. As stated earlier, this 
procedure implies that the NLO a{qq)‘^ correction provides the leading contribution in Ri . 

The various QCD parameters that will be used in the phenomenological analysis of (|^) are 

~ 0.3GeV , as{qqf = 1.8 x lO'^^GeV® , mq{qq) = (34) 

{asG'^) = 0.07GeV'‘ , /(G^) = 1.1 GeV^(asG^) , U = 0.132GeV , mq{qgGq) = l.5GeV'^mq{qq). (35) 

The parameter (agG^) represents the central value in the recent determination 0, 9^ (G^) is obtained from the 
dilute instanton gas model ||l8|| , {qgGq) is extracted from [l^ , and the dimension-six condensate parameter as{qq)^ 
is referenced to the vacuum saturation value which is known to underestimate the actual value by up to a factor of 2 
in the (/ = 1) vector and axial vector channels [ p0| . 

Before considering a detailed analysis of the sum-rules, we consider the Sq —> oo limit of the sum-rules which 
provides the following bound on the lightest resonance mass 


m2 < 


Ri (r, oo) 
Ro (r, oo) 


(36) 


which has the advantage of being a robust bound independent of the QGD continuum model. The explicit expressions 
for the sum-rules in the sq oo limit are 


Ro (r, oo) = 


1 


2407r^ TT 


a / 1301 a 

- 1 + 


1 


240 7r/r3 1207r2V7r/ 72 


/a\2 17 , 1 1 / 145 1 


1 16 a, , 47r / 1 a\ 47r /47 a > .o 

- - Y (1 + 108^ ) + Y ( 72^ ) 


i?i (r, oo) = 


1 


2407r^ TT 
1 

367r 


1 + 


1301 a\ 6 


1 - 


240 TT 
145 a 


1 


1207r2 




a\2 17 , ,1 

(11 - 67 b ) — 


72 


(aG^ 


1 16 a 
367r 9 TT 


(aG^ 


1 — 7b 


47r 

Y 


47 a 
72^ 


a{qqf 


(37) 


(38) 


The effect of the NLO {aG"^) and a{qq)'^ corrections is illustrated in Figure ||, where it is observed that the mass 
bound is increased significantly when these higher-order corrections are included. For brevity, we respectively refer 
to sum-rules containing the NLO and LO (aG^) and a{qq)'^ corrections as the NLO and LO sum-rules. In particular, 
we see from Fig. || that the /3(1600), excluded for the LO case, can be accommodated when the NLO corrections are 
included. Furthermore, the minimum of the NLO bounds occurs at a reasonable energy (r) scale in contrast to the 
rather large energy scale occurring when only LO corrections are included. 



FIG. 3: The ratio Ri (r, oo) /Ro (r, oo) as a function of r for the NLO (solid curve) and LO (dashed curve) sum-rules. As 
discussed in the text, acceptable values of Mr must lie below these curves. 

The mass estimate Mr is obtained by optimizing the choice of sq such that the most stable Ri/Rq ratio is obtained. 
Figures ^ and ^ illustrate this ratio for selected values of sq, resulting in Mr ss 1.6 GeV and sq « 4.0 GeV^ for the 
NLO case, while the LO analysis results in Mr ss 1.3 GeV and sq ~ 3.0 GeV^. These optimized values of Mr are 
consistent with the bounds established in Fig. ^ and explicitly demonstrate that the NLO condensate effects raise 
the estimated value of the hybrid mass. 
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FIG. 4: The NLO sum-rule ratio Ri (r, so) /Ro { t , so) as a function of r for the sequence of so values so = 

{3.0, 4.0, 5.0, 6.0} GeV^. The lowest (solid) curve corresponds to so = 3.0 GeV^, and the upper (dashed-dotted) curve 
corresponds to so = 6.0 GeV^. 



FIG. 5: The LO sum-rule ratio Ri (r, so) /Ro (t, so) as a function of r for the sequence of so values so = 

{2.0 ,3.0, 4.0, 5.0} GeV^. The lowest (solid) curve corresponds to so = 2.0 GeV^, and the upper (dashed-dotted) curve 
corresponds to so = 6.0 GeV^. 


The stability and self-consistency of the sum-rule analyses can be examined by using the optimized sq and Mu as 
input into the following expressions resulting from the single resonance model. 


.-Mir 


(t. So) — 
Ft 


(39) 


Figure || displays the left-hand side of ( ^^ for the NLO A: = 0,1, 2 sum-rules. We see from Fig. ^ that a stable ratio 
containing a r critical point occurs for each value of k, and that the variation of the corresponding value of JrM^ 
with k is minimal. By contrast, the corresponding curves for the LO sum-rules shown in Figure do not exhibit a 
critical point in the same r range associated with the Fig. mass estimate, and show strong dependence on k. Thus 
the NLO {olG"^) and a{qqY corrections lead to improved stability and self-consistency in the sum-rule analysis. 

The dominant uncertainties in the sum-rule analysis associated with the QCD parameters (34) and (^) are related 
to and a{qq)^ If the pattern established in other sum-rule channels is upheld for the hybrid OPE, then 
a{qq)‘^ given in ( p4|) underestimates the true value. Similarly, Aj^ = 300 MeV is a lower bound following from 
tts {Mr) iQ. Increasing either of these parameters increases the hybrid mass estimates, and thus it appears difficult 
to accommodate the /3(1400). 

Finally, we have verihed that the results of our analysis are essentially independent of the choice of renormalization 
scale = l/r motivated by renormalization-group improvement in the sq oo limit |p^ . Choosing renormalization 
scales in the energy region near the hybrid mass has a minimal effect on the sum-rule analysis. 


V. SUMMARY 

The (NLO) as corrections to the {agC^) and a{qq)‘^ contributions in the two point correlator of the current 
gq'yyiG^JT°‘q{x) have been calculated, and the effect of these contributions on the QCD sum-rule estimates of the 














0.507 



FIG. 6: The quantity Rk (r, so), which in the single-resonance model corresponds to is displayed as a 

function of r for the fc = 0, 1,2 NLO sum-rules. The optimized values so = 4.0 GeV^ and Mn = 1.64 GeV are used as inputs, 
the lowest (solid) curve corresponds to fc = 0, the intermediate (dotted) curve corresponds to fc = 1, and the upper (dashed) 
curve corresponds to fc = 2. 



FIG. 7: The quantity Rk (r, so), which in the single-resonance model corresponds to f%M%, is displayed as a 

function of r for the fc = 0, 1, 2 LO sum-rules. The optimized values so = 3.0 GeV^ and Mr = 1.31 GeV are used as inputs, the 
lowest (solid) curve corresponds to fc = 0, the intermediate (dotted) curve corresponds to fc = 1, and the upper (dashed) curve 
corresponds to fc = 2. 


1 ^ hybrid mass have been examined. The NLO a{qq)‘^ corrections are particularly interesting since they provide 

the leading contributions to the Ri sum-rule. The NLO (osG^) and a{qq)‘^ contributions improve the stability and 
self-consistency of the sum-rule analysis, resulting in a 1 ^ hybrid mass of approximately 1.6 GeV. This result reflects 
a lower bound devolving from the QCD input parameters, so it appears difficult to accommodate the /5(1400) as a 
hybrid state. 
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